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Abstract 

We perform an analysis of the relation between the factorization 
scale and the masses of the quarks in the calculation of the hard gluon 
coefficient in polarized deep inelastic scattering. Particular attention 
is paid to the role of strange and charm quarks at finite momentum 
transfer. It is found that for the momentum transfer of the present 
experiments, the contribution from the charm quark is significant. 

1 Introduction 

In the usual analysis of the proton spin structure function, based upon QCD 
and the operator product expansion (OPE), the moments of the singlet part 
of gi{x, Q"^) are written as: 



[\['\x,Q')x-'dx = AS4/.2)C^(QV/z^) + A^7„(/z2)C^(QV/i2), (1) 
Jo 

with C^i^Q"^/ fi^) = 1 + 0{as) the Wilson coefficients for the quark operators, 
C^{Q'^/fi'^) = 0{as) the Wilson coefficients for the gluon operators. The 
matrix elements AI]„(/x^) and A(7„(/i^) are not determined by perturbative 
QCD and should be either fixed by experimental constraints or calculated 
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using non-perturbative techniques. Eq. (|T]) can be inverted, using the inverse 
MeUin transformation, and the result is: 

where ® denotes a convolution of the two functions. 

Much of the debate on the proton spin in the last few years has been 
centered on whether or not the spin of the proton receives a contribution 
from the gluons [1-6]. On the basis of the OPE the picture is clear: there 
is no twist two gluon operator contributing to the first moment of gi, and 
hence Jq C^{x,Q'^ / n'^)dx = 0. This result implies that the first moment 
of Qi is given solely by the first moment of AE. If AS were identified 
with the spin in the proton carried by the quarks then the gluons would 
give no contribution. In this scenario, following the parton model language, 
AE„(/x^) = N J^f'^fnifJ''^), with A/ = A/i the amount of spin carried by 
the / quark and equals 1/9 for three flavors, 5/36 for four flavors, etc. It 
happens that AS cannot be identified with spin because of the axial anomaly. 
Indeed, the axial anomaly is at the heart of the disagreement between the 
OPE and the improved parton model (IPM) results for the role of gluons in 
the first moment of gi{x, Q^). In this contribution, we are not going to make 
a complete analysis of the equivalence (or otherwise) of these approaches but 
will limit ourselves to the analysis of the gluon contribution in the light of 
the IPM only. 

In the IPM the situation is more complicated. One calculates the full, 
polarized photon-proton cross section and uses the factorization theorem to 
separate the hard and soft parts: 

a^'''^{x,Q') = al''\x,Qyfi')^Af,/,,{x,f,') + al''<^{x,Qyfi')®Af,/^^^^ 

(3) 

where /x^ is the factorization scale, Afg(^gyj^ is the polarised quark (gluon) spin 
distribution inside the nucleon and at is the polarized, hard photon-quark 
or hard photon-gluon cross section. One then could relate gi calculated in 
the IPM, Eq. (^, to gi calculated in the OPE, Eq. (||), by identifying the 
hard, perturbatively calculated, Wilson coefficients with the hard photon- 
quark and hard photon-gluon cross sections and identifying the matrix el- 
ement AS(a;,/i^) {Ag{x, fi"^)) with the factorized quark (gluon) distribution 
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A/g(g)/7v. However, as already mentioned, AS(x,/i^) cannot be identified 
with the quark spin distribution. The relation between them is beyond the 
scope of the present work. Instead, we will concentrate on the relation be- 
tween the Wilson gluon coefficient and the hard gluon cross section of the 
IPM. Although there are excellent treatments of this subject in the literature 
[0, H, 0, |[, we think that the present contribution adds significantly to the 
understanding of the behaviour of gi{x, Q"^) at finite Q^. 

In section 2 we will develop the basis for the calculation of the hard gluon 
coefficient in the IPM. The resulting expression interpolates the known limits 
of —^Nf for << fi"^ and for >> yU^ and overcomes convergence 
problems found in an early work [^. The effects of this generalized, hard 
gluon coefficient are discussed in section 3. In particular, its effect on the 
contribution to gi{x, Q"^) from up, down, strange and charm quarks is studied. 
Our results indicate that these corrections are sizable and must therefore be 
taken into account when extracting the polarized gluon distribution from the 
proton. We also point out in this section how this anomalous contribution is 
affected by finite Q^. In section 4 we calculate the amount of polarized gluon 
in the proton necessary to explain the available data. We compare our result 
with other estimates made using simply the limiting cases for the hard gluon 
coefficients. Section 5 is used to study the region in x where this contribution 
is located. In section 6 we summarise the results obtained in this article. 



2 Theoretical Construction 

The hard gluon cross section is extracted from the full photon-gluon fusion 
cross section, a"'"^ and is calculated through the box graphs which start at 
order as- The other contribution from which it must be separated is the 
quark distribution inside the gluon 0. Mathematically this is expressed as: 

Q') = ^T'ix, QV/i') + Ag^^, f^'), (4) 

where Aq^ is the polarized quark distribution inside a gluon and ctJ[ ^ is the 
hard photon-gluon cross section defined, in the IPM, as the contribution 
coming from quarks in the box graph with transverse momenta greater than 
the factorization scale. 

The full photon-gluon cross section has been calculated to be M, ITD[ : 
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a 



-—Nf 



4^2 p2 



(2^-l)(l-^) 




with = — the gluon virtuahty, niq the quark mass and = ^ i'^-^)-P ^ 
the invariant mass squared of the photon-gluon system. For very large mo- 
mentum transfer, » 'm'^,P'^, the full cross section reduces to: 



—N 

2tx ■ 



(2x 



In 



+ P'^x{l — x) 



X] 



2ml - ^^^3^(20; - 1) 



+ P^a;(l — x) 



In- 



X 



(6) 



It remains to calculate Ag^. This is given by computing the triangle diagram 
or, equivalently, the integral over the transverse momentum of the square of 
the norm of the light-cone qq wave function of the gluon p, 0, |Tl]] . As Ag^ is 
a soft contribution, the integral over the transverse momentum has to have 
a cut off: 



Aq'^(x,U ) = —Nf / dfc, r p5-pr 



a 



-^Nf{{2x-l)l 



n 



'/i^ + P^x(l — x) + m^' 



+ (1 ^) 2m^ + PMl-2x) | ^ 

+ p2x(l — x) jj? + P2x(l — x) + J 

Equation is a generalization of previous results 0, including the depen- 
dence on the factorization scale for any values of the quark masses and gluon 
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virtuality. Its first moment is zero for /i^ << m'^,P'^. If /x^ >> m'^,P'^ the 
first moment of Aq^{x, jj?) is for P^ » m^, while it is ^Nf for >> P^. 
Using Eqs. (^), (^) and (|^ we can calculate the hard gluon coefficient: 




2ml + 



In ( — — ^ I + In ^ ^ 



^2 _|_ p2^(^i _ _|_ ^2 y \ X 



+ ^ 




^2 _j_ y^2 _|_ p22;(^]^ _ X 



Notice that the first moment of Eq. (P) does not depend on the ratio / P^ in 
the region /i^ >> m^, - it is a legitimate hard contribution. Equation (|^) 
is also a generalization of previous results and from its limit, /i^ >> m^, P^, it 
may be argued that the gluons contribute to the first moment of gi{x, Q^) 
because Jq crl ^{x,Q'^)dx = —^Nf. 

On the other hand, if one calculates the quark distribution inside a gluon 
through the triangle graph, which we call AqQp^, using a regularization 
scheme that respects the axial anomaly, it is found that Q: 




x] 



Ij2 _|_ p2x(^i — x) + m? 



(9) 



where the renormalization scale in the regularization of AqQp^ (using MS) 
has been taken to coincide with the factorization scale in the IPM. 

Equipped with Eq. (|^) we can understand exactly why the hard gluon 
coefficient in the IPM has a first moment different from zero. The reason 
is that in the process of factorization the axial anomaly was shifted from 
the quark distribution inside the gluon to the hard coefficient. Equation (|^) 
reflects the fact that the regularization of Aqpp^ respects the axial anomaly 
while the regularization of Aq^ does not. We also see that in the limit >> 
/i^, the discrepancy between the two calculations disappears (at least for the 
first moment - the x dependence depends on the regularization method). A 

^The triangle graph regularized with a cut off on the transverse momentum results in 
Eq. (|). 
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similar phenomenon is found in unpolarized deep inelastic scattering where 
an analysis by Bass |12] has shown that the trace anomaly induces the same 
sort of shift when a cut off over the transverse squared momenta of the quarks 
is used to separate the soft and hard regions. 

As a consistency check of our equations, we calculate the OPE hard co- 
efficient, . It is defined in the same way as aj^ ^ in Eq. (|) and calculated 
with the help of expressions (P) and (H): 



{2x-l)(ln% + ln^^ ^ 



X 



1 +2(1 -xl 



(10) 

This result is independent of mass and its first moment is always zero, in 
accordance with the results of Kodaira 13| and Bodwin and Qiu 0]. 



Consequences for the First Moment of the 
Hard Gluon Coefficient 



It is interesting to study the dependence on fi of the first moment of 
In an early study on this subject^, Mankiewicz and Schafer determined 



the first moment of the box graph as a function of the minimum transverse 
momentum carried by the quarks. Their results for Q"^ ^ oo agree quali- 
tatively with ours, as will soon be seen. But it was also found in Ref. 



that for momentum transfers of the order of 10 to 100 GeV^, the contribu- 
tion from light quark^ {rriq = 10 MeV) is deeply affected by the choice of 
the minimum value for the transverse quark momentum. In the method used 
here, such an ambiguity does not exist for the light quarks and its anomalous 
contribution for = 10 or 100 GeV^ is well defined and independent of k±. 
We use this result to argue that the hard gluon coefficient calculated here is 
more stable from the point of view of infrared singularities. 

Even with the known variations of the anomalous contribution with the 
factorization scale, it has been widely assumed in the literature |§ that for 

^We thank S. Bass for pointing out to us this work. 

■^We assume for the quark masses their current values. We do not take into account 
variation of the masses with the factorization scale but note that our conclusions are not 
significantly altered by small changes in the quark mass. 
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light quarks (u, d and s) the first moment of a1 ^ is —as/2n and for heavy 
quarks (hke c or b) it is zero (because, for >> /i^, ^ reduces to C^). 
But it also happens that the gluonic contribution to gi{x, Q"^) is of the form 
a1 ^{x, Q^/ ji^) ® A(7(x, jj?). This means that the scale jj? at which the gluon 
distribution is calculated (or parametrized) is the same scale /i^ that has to 
be used in the calculation of the hard gluon coefficient, and that it does not 
make sense to talk about the magnitude of the hard gluon coefficient without 
specifying the factorization scale. Thus, the heavy quark contribution is 
negligible only when the polarized gluon contribution is calculated at a very 
low scale compared with the quark mass. 

In Fig. m we show the first moment of function of the factor- 

ization scale for the u and d quarks (m^ ~ 25 x 10~^ GeV"^), for the s quark 
(m^ ~ 0.04 GeV"^) and for the c quark (m^ ~ 9/4 GeV"^). We see that, as is 
well known ^, 0; l^], the c quark does not contribute when ml » fi"^, as 
one can also verify directly from Eq. (^). However, for reasonable values of 
fi"^ there is a contribution large enough to be taken into account. Thus, the 
significance of the charm contribution to gi{x, Q"^) depends on where the po- 
larized gluon distribution is calculated. For instance, calculations have been 
made in the literature using input polarized gluon distributions at a scale 
of typically 4 GeV^. The authors of these calculations usually disregard the 
charm contribution. We note in passing that in the region of //^ where polar- 
ized charm can be disregarded, the polarized strange quarks yield only half of 
the contribution given by u and d quarks. As we see from Fig. |l|, the c quark 
gives around 64% of the contribution of the light quarks for /i^ ~ 4 GeV^ 
and so it should not be disregarded if the gluon distribution is calculated at 
this scale. We also see from Fig. ^ that the u and d quarks give the same 
contribution, independent of the factorization scale. We further notice that, 
for practical purposes, the hard gluon coefficient is independent of the exact 
value of the gluon virtuality P^. 

The discussion of the preceeding paragraph was based on the not so re- 
alistic assumption that the momentum transfer is infinitely bigger than 
any other scales in the theory. It implies, for instance, that when inte- 
grating the hard cross section we allow x to go from zero to one. But 
from simple kinematic arguments we know that x has a maximum value 
of Xjnax = Q^/ iQ"^ + + 4m2) and so Xmax 1 only when Q'^ » m^, P^. 
For the finite of the current experiments, Xmax never reaches 1 and so the 
integral in x has a cut off. For instance if one calculates the first moment 
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of o\ ^ for the c quark {m\ = 9/4 GeV'^) at /i^ = = 4 GeV^, using Eq. 
(H), one finds that its value changes from -0.64, when x is artificially allowed 
to reach 1, to ~ 0.015 when the physical cut off in x is applied. What hap- 
pens is that expression (|) itself was obtained under the assumption of an 
infinitely large . To be more consistent when dealing with finite Q^, one 
should derive the hard cross section from the full cross section without any 
approximation. 

In the general case we then write: 



C = cj^-' - ^q'oPE, (11) 
^ra ^ ^7.9 _ Ag9, (12) 

with a"^"^ given by Eq. and Ag^ and Aqpp^ given by Eqs. (|^ and (|^). 
We stress that these equations are the complete result at order a^. In Figs. || 
and ^ we show the first moment of ctJ[ ^, defined in Eq. (12), as a function of 
the factorization scale /x^ for = 10 GeV^ and = 3 GeV^, respectively. 
These values were chosen because they are the average of the EMC [|19|, ^ 
and SLAC |^ experiments. The resulting dependence is very interesting. It 
shows that in the region of interest (yU^ > 1 GeV"^) there is no appreciable 
dependence on the gluon virtuality or on /^^ (at least for = 10 GeV^) 
but the mass dependence is strong. Remarkably, the contribution from the s 
quark is never the same as the contribution from the u and d quarks, contrary 
to what is usually claimed. The s contribution is ~ 0.9 as{Q^) /2tt for the 
EMC data and ~ 0.75 as{Q^)/27i for the E143-SLAC data. We also find 
a nonnegligible contribution coming from the c quarks. For the EMC data, 
the c quark contributes with ~ 0.2 as{Q'^)/'2'K and for the E143-SLAC data 
with ~ 0.1 as((5^)/27r. Figure ^ is unaltered ^ if we go from = to 
rrig = 1 X lO^^Ge^^. If we then compare our Fig. |^ with Fig. 2 of Ref. 

we see clearly that the present approach does not have a convergence 
problem in and yelds a perfectly unambiguous contribution from the light 
quarks. Finally, we show in Fig. ^ the dependence of the polarized charm 
contribution calculated with /x^ = 3 GeV^. This contribution, obviously, 
tends to the value calculated in Fig. [I| (~ —0.57 asiQ"^) /2-n). 

^In reality, there is a ~ 1% correction for /i^ ^ 0.01 GeV^ . This result is in complete 
accord with the fact that the anomalous contribution for the u and d quarks goes to zero 
as ix^ goes to zero. 
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4 Relevance in Analysing the Fraction of Nu- 
cleon Spin Carried by Gluons 

In terms of the polarized quark and gluon distributions, gi{x, Q"^) for 4 flavors 
is written as: 

+ 4AS(a;, + ^crl"^ix, Q^ /x^) ® Agix, /i^), (13) 
3d 3d 

where Agg = Au- Ad, Aqs = Au + Ad- 2 As, Aqi^ = Au + Ad + As- 3 Ac 
and AS = Am + Ad + As + Ac. For 3 flavors the coefficient of the singlet 
part changes from 5/36 to 1/9 and Agis does not exist. To order aslQ"^) 
the first moment of ( |T3|) is: 



nm = hm + hm - h,m + w 




2^ 



...^ + c,-M)^„g.). (14) 

Zn Zn j 

The coefficients of a<j(Q^) have the following meaning. The 2 indicates that 
the u and d quarks give the same contribution ""^^ \ as discussed before. The 
S\ and C\ factors give the amount of strange and charm quark contributions, 
according to Eq. ( [T^ ) and Figs. 1-3. 

To extract the value of AG(Q^) we will closely follow Refs. [0, [l^. For 
the sake of comparision, we begin with only 3 flavors and with the common 
assumption that the u, d and s quarks give the same anomalous contribution. 

Under the assumption that the polarized sea originates exclusively from 
the anomalous gluon contribution we have, for 3 flavors, the following iden- 
tities: 



h 

h + Io 



12 
1 

36 



{F + D){1 
(3F - D) 



vr 



TT 



+ 4(1 - 



37r' 



(15) 
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where the the quark spin fractions were expressed in terms of the F and D 
couphngs and corrections from the two loop expansion of the beta function 
and anomalous dimension were incorporated. In NLO, ag is given as the 
solution of the following transcendental equation: 



In 



A2 poa. 



A, 

^0 



Air 



+ 



(16) 



with /3o = 11 -2A^//3 and /?i = 102-38iV//3. We use A = A^^) = 248 MeV, 
determined by fixing A'-'^-* = 200 MeV |T8[. Using the experimental values 
of F and D as given in [0 , we determine I3 and Jg + Iq at = 10 GeV'^ 
(with a,(g2 = 10 GeV^) - 0.209): 



/3 = 0.0977 ±0.001 
Is + Io = 0.0779 ±0.002 (17) 

We now use Eq. (|l^) to determine AGIQ"^). On the left hand side, we use 
the experimental result P(]| : 

rP(Q2 = lOGeV^) = 0.142 ± 0.008 ± 0.011. (18) 

On the right hand side we use the results (|l3), Si = 1, Ci = and remember 
that for 3 flavors the singlet coefficient is 1/9 and /15 = 0. The result is: 

A^(g2 = lOGeV^) = 3.04 ± 1.4. (19) 
For 4 flavors the analysis is similar. One just has to redefine the integral of 

We then proceed as before and calculate AG using the result for the 
gluon coefficient as displayed in FigJ^. We see that for jj? = Q"^ = 10 GeV^ , 
si ~ 0.9, ci ~ 0.21 and a,(g2 = 10 GeV^) = 0.2142, resulting in: 
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h = 0.0976 ±0.001 
/8 + /0-/15 = 0.0786 ±0.002 
A^(Q2 = 10 Gel/^) = 2.32 ±1.06 (21) 

In passing we notice that if the usual assumption of infinite momentum trans- 
fer were used, then according to the results of Fig.[l|, at 10 GeV^ Si = 1, 
Ci ~ 0.81 and hence Ag ^ 1.89. 



5 The X dependence 

The exact x dependence of the anomalous contribution is a matter of con- 
vention because the freedom in the factorization scheme while calculating 
cl"^. Other choices of regularization would result in different functions of 



X. But, as shown by Gliick et al. |]14|, the exact form of the x dependence 



seems not to be very important. Once we do not know the form of the po- 
larized gluon distribution, the best we can do is constrain it by some general 
considerations. For instance, there is the positivity condition: 

\Ag{x,Q')\<g{x,Q^), (22) 

where g{x, Q"^) is the unpolarized gluon distribution. A very simple form 
that satisfies the above condition is: 

A(x)=x°^(x), (23) 

where a is determined through the normalization of Ag. For Ag of Eq. (|2T|) , 
a = 0.49. The advantage of using this form to study the x dependence is its 
simplicity. The problem with Eq. (p3D is that it does not have the correct 
behavior a.s x 0. As proposed by Brodsky et al. 



^ - (24) 

g{x) 

as X — > 0. From the many ways to satisfy both conditions (^) and (|2^) , we 
choose: 

Ag{x, /x^ = 9 GeV'^) = axg{x, jj^ = 9 GeV^){l - xf, (25) 
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where a = 6.92 for Ag = 2.32. We made this choice guided only by the desire 
of simphcity and to produce a polarized gluon distribution that resembles an 
already existing one |T5[ . For the unpolarized gluon distribution, we use the 



one given by the NMC p4|, determined from inelastic J/ip production: 



xg{x) = ^{r] + l){l-xr. (26) 

This parametrization is valid for /i^ = 9 GeV^ and should not be trusted for 
X < 0.01. Again, this choice is based on simplicity and we note that a further 
parametrization by the NMC group agrees with Eq. (^Bj) for x > 0.01. 
The parameter rj is rj = 5.1 ± 0.9. Given these choices, we show in Fig. ^ 
the forms (0) and (|25| ) for the polarized gluon distributions plus the forms 



of Brodsky et al. [|23[ and Gehrmann and Stirling (GS) |T^, calculated at 
4 GeV^. Our parametrization (|25|) is slightly higher than that of GS because 
of the normalization factor. If we use the same normalization as theirs^ both 
curves would be essentially the same. Evolution from 4 to 9 GeV^ for the 
GS distribution also has small effects. The parametrization of Brodsky et al. 



2^ is much smaller than the others because in their approach the polarized 
gluons are not responsible for the small experimental value of Eq. (p^Sf). 

Using the constructed gluon distributions, we can estimate where in x 
the anomalous contribution is located. In Fig. ^ we show the anomalous 
contribution, ^cr^ ^ {x , , fj,"^) ® Ag {x , fj,"^) , for Ag{x) = axg{x){l — x)^. We 
see that its contribution inside the experimental region is important. To better 
evaluate its importance, we calculated the amount of the total gluon that lies 
inside the region 0.01 < x < 1. For 4 flavors, the contribution from x > 0.01 
corresponds to about 66% of the total anomalous contribution. In the case 
of 3 flavors, this percentage is ~ 69%. For Ag{x) = x"g{x) this conclusion 
is not dramatically altered. 

It is also interesting to compare the anomalous gluon contribution directly 
with the experimental data. To this end, we plot in Fig. |^ the experimental 
data |1^, pO[ for gi{x) together with an early next-to-leading order estimate 



2^ for the valence quark distribution and the anomalous gluon contribution 
for the case of 3 and 4 flavors. A remark is necessary here. The calculation of 
the hard gluon coefficient was performed through a cut off on the transverse 



We note that in [y_5| , the couphng constant is calculated in leading order rather than 
in NLO. This would lead to an increase of the total polarization carried by the gluons. 
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momenta of the partons in order to regularize the integrals. This procedure 
is the definition of the parton model. On the other hand, we calculated the 
strong coupling constant, and also the evolution of g{ in Ref. ||2^, using 



the MS scheme. In principle, showing the x dependence of two quantities 
in different schemes is not a consistent procedure. However, the problem 
is not as bad as it looks. First, if we change schemes we can mantain 
unaltered by a simple redefinition of the parameter A in expression ([T6|) . 
Second, the theoretical curve for gi{x, Q"^) is to be interpreted as a guide 
of what a parametrization for the valence part of the polarized structure 
function would give, the regions in which it differs from the data and where 
it should be corrected. A proper procedure would be to calculate the quark 
distribution, the anomalous dimensions and the Wilson coefficients in the 
same scheme. That said we proceed noticing that the integral over x of the 
valence contribution calculated in (~ 0.169) is in complete agreement 
with the estimates calculated previously in Section 4. The two curves below 
the origin, are the anomalous contributions that should be added to the solid 
curve for Nf =3 or 4. As we fixed the normalization of the total polarized 
glue for either 3 or 4 fiavors, there is no noticeable difference between the 
two cases. We see that the anomalous contribution is potentially important 
to correct the x dependence of the polarized valence distribution inside the 
proton. 



6 Discussion 

In summary, there is a gluonic contribution to the proton spin when the 
IPM hard gluon coefficient is defined through Eq. (JT^) with Ag^ defined as 
the quark distribution inside the gluon with transverse squared momentum 
less than the factorization scale. As a consequence, this anomalous gluonic 
contribution in the IPM is free of infrared ambiguities. We showed that if 
we accept the commonly used assumption of an infinitely big momentum 
transfer, there is a c quark contribution to the spin in addition to the u, 
d and s quark contributions. The contribution from the massive quarks is 
dependent on the factorization scale at which the polarized gluon distribution 
is calculated. The c quark contribution is small only in the region /i^ < 
1 GeV^, in which case the s quark contribution is also strongly affected. 
We also calculated what would be the possible anomalous corrections 
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when the momentum transfer is in the region of the present experiments. To 
perform such a calculation we have to keep all terms in rn^/Q^ and P^/Q^ 
in the full photon-gluon cross section when calculating the hard gluon coef- 
ficient. This means that we are including higher twist effects and, although 
we use the complete result at order ^^(Q^), possible corrections coming from 
higher order terms in ^^(Q^) could be important and so our calculation is in- 
complete. Even so, we think that our results are more consistent than simply 
using the approximate expression for the hard gluon coefficient in the case 
of massive quarks and relatively low Q^. The corrections due to finite are 
not small and we think they should be taken into account when calculating 
the amount of spin carried by gluons. When studying the x dependence of 
the anomalous contribution, we conclude that both 3 and 4 flavors give ap- 
proximately the same contribution inside the experimental region. But the 
amount of polarized glue needed to fit the data is much smaller when charm 
is included. Moreover, we showed that from the conceptual point of view, it 
would be wrong not to include a fourth flavor. 

We would like to thank to S. Bass, S. J. Brodsky, W. Melnitchouk and 
G. Filler for helpful discussions. This work was supported by the Australian 
Research Council and by CAPES (Brazil). 
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Figure 1: Hard gluon coefficient as given by Eq. (§), calculated with the 
assumption of infinite momentum transfer as a function of the factorization 
scale. For realistic scales (/x^ > 1 GeV"^), the charm contribution is seen to 
be important. 
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Figure 2: Hard gluon coefficient as given by Eqs. (|), (|^ and (|T2|). The 
momentum transfer is fixed at 10 GeV^. It is seen that the strange quark 
contribution never equals that from the up and down quarks and the charm 
quark contribution is sizable for /i^ > 1 GeV^. 
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Figure 3: Hard gluon coefficient as given by Eqs. (|), (|^ and (|12D. The 
momentum transfer is fixed at 3 GeV^. It is seen tliat tlie strange quark 
contributes witli approximately 75% of the up and down quarks in the re- 
ahstic region of /x^ > 1 GeV^. In the same region, the charm quark gives a 
10% contribution. 
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Figure 4: Hard gluon coefficient for the charm quark, calculated with Eqs. 
(§)' (0) (0)- "^^^ factorization scale is fixed at 3 GeV^ and the 
dependence is studied. 
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Figure 5: Comparison of various polarized gluon distributions considered in 
the text. 



21 



— 010— contribution 

_ _ _ _ Up + down contribution 

Strange contribution 

Charm contribution 



-0.01 5 I , , , , , II , , , , , 

0.01 0.10 1.00 

X 



Figure 6: Comparison of the x-dcpendence of the non-strange, strange and 
charm quark distributions to gi{x). The anomalous contribution is given by 
^{x, Q'^, iJ?) ® Ag{x, /x^) and it is used the form Ag{x, /x^ = 9 GeV'^) = 
axg{x){l — x)^ for the polarized gluon. 
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Figure 7: EMC |T9[ and SMC 13 data for gip{x) at 10 GeV\ The theoretical 
curve for the polarized valence distribution is calculated in NLO and taken 
from Ref . ||26[ - 



The anomalous contribution should be subtracted is to be 



subtracted from the theoretical curve. 
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